第 38 卷 第 3 期 工程 热 物 理学 报 
2017 年 3 月 JOURNAL OF ENGINEERING THERMOPHYSICS Mar., 2017 


AN 人 全 匡 日 
LU FY 


Vol.38, No.3 


相 变 换 热 系统 的 能 量 流 模型 及 其 应 用 
陈 三 1 齐 腾 云 2 着 春 涛 3 钟 英 席 3 起因 玉 ! 陈 群 


(1. 清华 大 学 工程 力学 系 , 热 科学 与 动力 工程 教育 部 重点 实验 室 ， 北 京 100084; 
2. 中 电 国际 新 能 源 控股 有 限 公 司 ， 上 海 “200086; 
3. 东莞 中 电 新 能 源 热电 有 限 公 司 ， 东 莞 523127) 


摘 要 工 质 在 换 热 过 程 中 发 生 相 变 时 , 其 物性 参数 将 发 生 剧 烈 变 化 , 增 大 了 换 热 系统 性 能 分 析 和 优化 的 复杂 度 . 本文 以 
顺 流 、 道 流 和 了 叉 流 三 种 典型 换 热 避 为 研究 对 象 , 在 考虑 工 质 发 生 相 变 的 前 提 下 , 构建 了 与 其 对 应 的 等 效 热 阻 网 络 图 (能 量 
流 模 型 )， 获 得 了 换 热 器 的 换 热 量 与 结构 /运行 参数 以 及 工 质 物性 间 的 直接 物理 联系 . 以 此 为 基础 ， 对 于 含 相 变 工 质 的 换 热 
系统 ,通过 构建 其 整体 能 量 流 模型 ， 结 合 基 尔 霍 夫 电压 、 电 流 定律 明确 了 系统 中 各 独立 设计 参数 间 的 整体 约束 关系 ,并 应 
用 拉 格 朗 日 乘 子 法 实现 了 相 变 换 热 系统 的 性 能 优化 . 


关键 词 换 热 系统 ; 相 变 ; 能 量 流 模 型 爆 耗 散热 阻 ; 性 能 优化 
中 图 分 类 号 : TK121 文献 标识 码 : A ”文章 编号 ;0253-231X(2017)03-0600-07 


Energy Flow Model for Heat Transfer Systems With Phase Change Fluids and 


Its Applications 


CHEN Xil QI Teng-Yun? JIANG Chun-Tao? ZHONG Ying-Biao? ZHAO Si-Yu! 


(1. Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, 
Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China; 
2. China Power International New Energy Holding Ltd, Shanghai 200086, China; 

3. China Power New Energy Dongguan Cogeneration Co., Ltd, Dongguan 523127, China) 
Abstract Phase change of working fluids during heat transfer processes leads to drastic variation of 
fluid properties, which makes the heat transfer analysis and optimization much harder. This study 
first introduced the equivalent thermal circuit diagrams, i.e. energy flow models, and deduced the 
mathematical relations among the design requirements, structural and operating parameters and fluid 
properties for parallel-flow, counter-flow and cross-flow heat exchangers with phase change fluids. On 
this basis, we constructed the energy flow model for a multi-loop heat exchanger network (HEN) with 
phase change fluids and then clarify the integral relations among each independent design parameter 
of the HEN with the aid of Kirchhoff’s voltage and current laws. Finally, by applying the Lagrange 
multiplier method, the performance of the multi-loop HEN was analyzed and optimized to show the 
applications of the newly proposed energy flow model. 
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0 引 各 


换 热 讨 程 广泛 存在 于 能 源 、 动力、 石油 、 化 工 等 
工业 领域 ， 其 性 能 的 优 劣 将 直接 影响 到 系统 的 整体 
性 能 ， 因 此 对 换 热 过 程 及 其 系统 进行 性 能 分 析 和 优 
化 设计 得 到 了 广大 学 者 的 长 期 关注 . 其 中 , 对 数 平均 
温差 法 和 e-NTU 法 等 经 典 方法 在 换 热 需 性 能 设计 和 
校 核 中 发 挥 着 重要 作用 .近年 来 , Chen 等 ,3 将 换 
热带 中 换 热 过 程 的 能 量 守 人 恒 方 程 和 传 热 方程 代入 燃 
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耗 散热 阻 的 定义 式 问 , 在 不 考虑 工 质 物性 变化 的 前 
提 下 ， 导 出 了 顺 流 、 道 流 和 又 流 等 典型 换 热 句 的 烽 
耗 散热 阻 的 计算 式 ， 并 给 出 了 单个 换 热 絮 的 等 效 热 
阻 图 . 同时 ， 针 对 多 回路 、 串 联 和 并 联 等 换 热 系统 ， 
构建 了 相应 的 等 效 热 阻 网 络 图 向 ， 提 出 了 反映 换 热 
系统 中 热量 整体 输 运 规律 的 能 量 流 模型 ， 并 结合 

尔 霍 夫 电压 、 电 流 定律 导出 了 换 热 系统 性 能 与 其 结 
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构 / 运 行 参数 以 及 工 质 物 性 间 的 整体 约束 关系 . 
考虑 到 换 热 过 程 中 工 质 的 物性 变化 , Zhao 等 5 
ss 
流 模型 ， 并 在 不 考虑 工 质 压 力 变化 的 前 担 下 ， 对 多 
回路 换 热 右 网 络 进行 了 性 能 优化 。 然而， 当 换 热 恬 
内 的 工 质 发 生 相 变 时 ,不 仅 工 质 物性 和 换 热 右 换 热 
系数 发 生 突变 外， 而 且 也 会 对 闭 式 换 热 系统 的 工作 
压力 产生 影响 。 因此， 本 文 在 考虑 工 质 相 变 和 压力 
变化 的 前 提 下 ， 分 别 对 顺 流 、 逆 流 和 又 流 换 热 需 建 
立 等 效 热 阻 网 络 模型 。 以 此 为 基础 ， 构 建 多 回路 换 
热 系 统 的 整体 能 基 流 模型 ， 结 合 基 尔 霍 夫 电 压 和 电 
流 定律 推导 工 质 压 力 变 化 时 ， 换 热 系统 中 各 独立 参 
数 间 的 整体 约束 关系 ， 并 应 用 拉 格 衣 日 乘 子 法 ， 对 

相 变 换 热 系统 的 整体 性 能 进行 优化 . 


1 相 变换 热 种 的 热 阻 网 络 模型 


1.1 含 相 变 过 程 的 顺 流 换 热 器 的 热 阻 网 络 模型 

图 1 为 热流 体 疡 与 冷 流体 ! 在 顺 流 换 热 右 中 换 
热 的 工 一 g 图 , 热流 体 进 、 出 口 温 度 分 别 为 7h,1 和 
Th,2, 冷 流体 进 、 出 口 温度 分 别 为 ,1 和 TI,2， 换 热 
需 的 换 热 面积 为 4. 冷 流体 在 换 热 带 中 吸收 热量 发 
生 液 -气相 变 , 相 变温 度 为 了 ie. 图 1 反映 了 换 热 需 
中 冷 、 热 流体 温度 和 换 热 量 之 间 的 关系 , 其 曲线 斜率 
为 对 应 流体 热 容 量 流 的 倒数 . 可 见 ， 在 换 热 过 程 中 ， 
由 于 低温 工 质 发 生 相 变 ， 其 比 定 压 热 容 在 相 变 前 后 
发 生 突变 ， 在 预 热 段 为 工 质 的 液态 比 定 压 热 容 ， 在 
过 热 段 为 工 质 的 气态 比 定 压 热 容 ， 在 相 变 段 工 质 的 
比 定 压 热 容 可 认为 是 无 穷 大 。 同样 ， 工 质 在 不 同 相 
态 下 对 应 不 同 的 换 热 系数 ， 

如 图 1 中 虚线 所 示 , 将 换 热 面积 均匀 划分 为 N 
段 , 并 按 沿 低温 流体 的 流动 方向 顺序 编号 为 1~ NN 
每 个 换 热 单元 的 面积 为 A/N， 换 热量 为 Qi。Thm, 
和 了 insz 分 别 为 第 ;个 换 热 单元 中 热 、 冷 流体 的 特征 
温度 , 即 该 单元 热 、 冷 流体 各 上 自 进 、 出 口 温度 的 算术 
平均 值 . 若 分 段 数 足 够 大 ， 则 分 段 后 每 个 换 热 单元 
中 流体 的 相 态 和 比 定 压 热 容 可 认为 是 定 值 。 对 于 每 
个 换 热 单元 , 燃 耗 散热 姐 与 换 热 量 和 算术 平均 温差 
的 关系 ， 及 燃 耗 散热 阻 与 换 热 面积 、 传 热 系数 和 流 
体 热 容量 流 之 间 的 关系 为 外 
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式 中 , 下 标 h 和 1 分别 表示 热 、 冷 流体 ，m 表示 算 
数 平均 值 ，i 表示 第 i 个 换 热 单元 ，Rs 为 烽 耗 散热 
阻 , 大 为 换 热 系数 ，4; = A/N 为 换 热 面积 ，G 为 热 
容量 流 ，G=mcp,m 为 质量 流量 ,cj 为 工 质 的 比 定 
压 热 容 . 


图 1 低温 工 质 发 生 相 变 的 顺 流 换 热 器 人 工 一 4 图 
Fig. 1 The 7T — gq diagram of the parallel-flow heat exchanger 
with phase change of cold fluid 


根据 方程 (1) 并 应 用 热电 比拟 方法 , 图 2 给 出 了 
每 个 换 热 单元 的 等 效 热 阻 图 . 可 见 , 冷 热流 体 间 的 换 
热 过 程 可 以 等 效 为 热量 Qi 在 热 势 差 (Thm,i 一 Timi) 
的 作用 下 流 过 热 阻 Rs, 的 过 程 ， 


2 
一 一 > 
Th Tm 
Re; 


图 2 单一 换 热 单 元 的 等 效 热 阻 图 
Fig. 2 The equivalent thermal circuit for a segment in the 
heat exchanger 


根据 能 晤 守恒 方 程 ,， 流体 在 第 i 个 换 热 单元 的 
算术 平均 温度 与 其 进口 温度 、 各 有 段 的 换 热量 和 热 容 
量 流 间 满足 如 下 关系 


Th,m,i = Th,1 — (总 


Tm -= Ti1 十 区 


! 吕 名 | (2) 


> | (G) 
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整个 换 热 过 程 中 , 流体 的 算术 平均 温度 为 


全 了 1 : 
Th bt A 多 
二 h,s 
硼 和 站 钙 1 由 0 
1,1 起 党 本 
Tim = 二 了 1 十 一 5 
1 2 (5) 


联 立 方程 (1)~(5), 可 以 得 到 整个 换 热 过 程 中 热 、 冷 
流体 的 算术 平均 温度 Th,m 和 Tm 之 间 的 数学 关系 


Th,m J Tim 二 一 Ep 十 (0 (Rn,i Li 十 局 一 E1 (6) 


其 中 时 
_1 Qi 
Ee 2 Grn Wm 
= i—1 
四 上 (OF 
Mi 2Gh QiGn,j (8) 
1 QO 
Ri; = 十 9 
, 2G1i 2 Qi (9) 
N 
1 i 
二 二 写 2 di (10) 


根据 方程 (6) 并 应 用 热电 比拟 法 , 图 3 给 出 了 
顺 流 相 变 换 热带 的 等 效 热 阻 网 络 图 。 


图 3 顺 流 相 变换 热 器 的 等 效 热 阻 网 络 图 [5 
Fig. 3 The equivalent thermal circuit for the phase-change 
parallel-How heat exchanger 


式 中 ,en、a1 为 附加 热 势 ， 其 中 ，sh 将 热量 的 热 势 
由 高 温 流 体 的 算术 平均 温度 Thm 提高 至 入 口 温度 
人 1 sl 将 热量 的 热 势 由 低温 流体 的 入 口 温度 ,1 提 
高 至 算术 平均 温度 Thm。Rbi、Rii 为 附加 热 阻 ,其 
中 ,Rn 使 热流 Qi 的 热 热 由 Th 1 降 至 Th mi, Ri 使 
1.2 含 相 变 过 程 的 逆流 换 热 器 的 热 阻 网 络 模型 

对 于 含 相 变 过 程 的 逆流 换 热 器 ， 同 样 可 以 将 换 
热 面积 按 相 变 工 质 流向 均匀 划分 成 N 段 , 并 通过 传 


热 方程 和 能 量 守恒 方程 ， 得 到 整个 换 热 过 程 中 流体 
的 算术 平均 温度 Th 和 TY xm 之 间 的 数学 关系 


Th,m es Tm 二 一 Ep 十 (OF (Rpi 四 Rigs 4 Ri1ii) CY (11) 


其 中 
1 1 kiA 1 1 | 
Ghpi Qi SPUN Ci Os 


fsi=— 3 RBA/1l 
A 


式 中 ,附加 热 势 sa、sl 和 附加 热 阻 R;、Rii 的 表达 
式 见 文献 回 .方程 (6) 和 (11) 具有 同样 的 形式 ， 所 
以 根据 方程 (11) 可 以 得 到 与 图 3 类 似 的 等 效 热 阻 网 
络 图 , 但 其 中 附加 热 阻 加 ,; 和 已,; 的 表达 式 不 同 . 
1.3 含 相 变 过 程 的 又 流 换 热 器 的 热 阻 网 络 模型 

对 于 冷 、 热 流体 均 不 混合 的 叉 流 换 热 顺 ， 将 换 
热 需 的 换 热 面 积 沿 冷 、 热 流体 的 流动 方向 分 别 划分 
为 N 份 , 相应 的 质量 流量 也 分 为 N 股 . 此 时 , 该 又 
流 换 热 器 等 效 为 N? 个 冷 、 热 流体 各 自 混合 的 又 流 
换 热 单元 的 组 合 . 与 顺 流 和 逆流 换 热 需 类 似 ， 可 导 
出 整个 换 热 过 程 中 热 、 冷 流体 的 算术 平均 温度 Ti 
和 Tw 之 间 的 数学 关系 


Thm—Tlm=—én+Qii (Rnijt+ Regijt+ Ri,j)—er (13) 


其 中 


(12) 


kiA 
1 eMemii 1 
Rs,ij = Be 
;2)] ev Shii 一 1 
kiA 
1 emToai4+1 N? (14) 
2G1ij om 一 1 人 


式 中 ,各 热 势 、 附 加 热 阻 的 表达 式 及 其 等 效 热 阻 网 
络 图 见 文献 [5]。 
2 相 变 换 热 系统 的 能 量 流 模 型 

在 第 1 节 中 ,以 三 种 典型 相 变 换 热 需 为 研究 对 
象 ， 建 立 了 相应 的 等 效 热 阻 网 络 模型 . 对 于 由 大 量 
换 热 右 构 成 的 换 热 系统 ， 可 以 以 已 构建 的 三 种 典型 
换 热 器 的 等 效 热 阻 网 络 图 为 基础 ， 按 照 换 热 器 间 的 
连接 关系 得 到 对 应 的 能 量 流 模型 . 

图 4 给 出 了 一 个 典型 的 多 回路 换 热 系统 的 结构 
示意 图 , 包括 换 热 器 1 和 2、 储 液 铅 、 泵 以 及 相应 管 
路 . 在 换 热 器 1 中 ,热量 由 高 温 流 体 释放 给 中 间 回 
路 流体 , 中 间 回 路 流体 吸收 热量 由 液态 变 为 气态 . 在 
换 热 器 2 中 , 热量 由 中 间 回 路 流体 释放 给 低温 流体 ， 
中 间 回 路 流体 由 气态 变 为 液态 . 在 循环 泵 的 上 游 设 
置 储 液 饶 ， 且 底部 安装 高 度 与 泵 和 换 热 器 的 高 度 一 
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致 . 工 质 由 储 液 饶 底 部 抽出 , 流入 循环 泵 , 使 循环 泵 
中 的 工 质 维持 在 液态 , 从 而 保证 循环 泵 的 安全 运行 . 
图 5 给 出 了 图 4 所 示 的 换 热 器 网 络 的 等 效 热 阻 网 络 
图 (能 基 流 模型 )， 其 中 下 标 1 和 2 分 别 用 于 表示 换 
热 恬 1 和 换 热 器 2. 


广 | EN 
图 4 多 回路 换 热 系 统 的 结构 示意 图 


Fig. 4 The sketch of a multi-loop heat transfer System 


二 wa FF ee | Ts 
、 a o | = 
和， 
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图 5 多 回路 换 热 系统 的 能 量 流 模型 [5 
Fig. 5 The equivalent energy flow model for the multi-loop 
heat transfer System 


基于 换 热 系统 的 能 流量 模型 ,热量 在 整个 换 热 
系统 中 的 输 运 规 律 就 可 以 由 炮 耗 散热 阻 以 及 基 尔 霍 
夫 “电压 ”和 “电流 ”定律 来 定量 描述 . 对 于 图 5 所 
示 的 能 基 流 模型 ，3 个 市 点 存在 2 个 独立 的 电流 方 
程 (15)、(16)，(2N+1) 条 文 路 存在 (2N 一 1) 个 独立 
的 电压 方程 (17)~(19). 


N 
>》， Q1i 一 Qtotal =0 


(15) 
4 
N 
> 2i Qtotal 一 0 (16) 
| 
Q1s (Ripit Rigi + Ri,ai) = 
Quit1 (Ripiti + Rigitl t+ Riait1) 
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Q2i (Ro2,ai + Ro,gi + R21i) = 


Q2it1 (Rozaitl tt Rogitl +t R21it1) 


(i=1,2,...,N-1) (18) 
Th,m 一 了 1m 十 El,h 未 E1,a 十 E2,a 十 E91 
Q11 (Rini + Ryg + Rial)— (19) 
Q21 (Roalt Rogl+t R211)=0 
其 中 
N 
enh = (QLk/Gunk)/2 (20) 
jl 
N 
ELa = (Qix/Guar)/2 (21) 
k=1 
N 
Eg = 》 (Q@2k/G2 ak)/2 (22) 
el 
N 
E21 = > (Qok/G218)/2 (23) 
k=1 
i 
Rini = 1/(2G1,ni) 十 >》， (Qk/ (QiGLn,k)) (24) 
B= 
N 
Riai=1/(2G1ai)+ > (Qin/(QiiGuan)) (25) 
天 一 ?十 1 
| 
Raai 王 1/(2G2 si) 十 (@2k/(Q2iGa ap)) (26) 
Xi 


N 
R21i = 1/(2G2,1i) 十 和 (Qo /(Q2iG2,1x)) (27) 


Ri,gi = (1/(2G1,ni) — 1/(2G1,ai)) x 

(1+2/(exp(kiiAi(1/Gini—1/G1,ai)/N)—1)) (28) 
Ro,gi = (1/(2G2,3i) — 1/(2G2,1i)) xX 

(1+2/(exp(k2,iA2(1/G2,a0) /G2,a,1/G211)1)) 


由 于 附加 热 势 sin 将 热量 的 热 势 由 高 温 流 体 的 算术 
平均 温度 提高 至 入 口 温度 ， 而 附加 热 势 = 将 热量 
的 热 势 由 低温 流体 的 入 口 温 度 提高 至 算术 平均 温度 ， 
因此 式 (19) 等 价 于 

hl1 一 全 1 十 sla 十 E2a 一 


Q11 (Rini +t Rug + Rial)— 
Q21 (Roalt Rogl+t R211)=0 


(29) 


(30) 


对 于 换 热量 给 定 的 换 热 器 网 络 ， 降 低 其 制造 和 
运行 成 本 通常 是 系统 优化 的 主要 目标 . 换 热 避 网 络 
的 总 换 热 面积 是 主要 制造 成 本 ,而 各 文 路 内 工 质 的 
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总 质量 流量 是 主要 运行 成 本 ， 因 此 对 于 换 热 器 网 络 
的 优化 目标 可 以 简化 为 换 热 器 的 总 面积 最 小 和 工 质 
的 总 质量 流量 最 小 . 通常 情况 下 ， 降 低 其 中 一 种 成 
0 种 成 本 的 增加 ， 因 此 换 热 器 网 络 的 
结构 和 运行 参数 的 优化 可 以 转变 为 以 下 两 种 优化 问 
题 : 1) 在 各 支 路 工 质 的 总 质量 流量 给 定 的 前 提 下 , 优 
化 换 热 面积 和 流体 流量 分 配 ， 使 得 所 有 换 热 句 的 总 
换 热 面积 最 小 ; 2) 在 所 有 换 热 器 的 总 换 热 面 积 给 定 
的 前 提 下 ， 优 化 流体 流量 和 换 热 面积 分 配 ， 使 得 各 
文 路 流体 的 总 质量 流量 最 小 . 对 于 图 4 所 示 的 换 热 
器 网 络 , 若 冷 、 热 流体 的 入 口 温度 及 比 定 压 热 容 , 换 
热量 ， 以 及 换 热 器 1 和 2 的 总 换 热 面积 或 换 热 工 质 
的 总 质量 流量 给 定 ， 通 过 优化 各 换 热 器 的 换 热 面积 
和 各 文 路 工 质 的 质量 流量 可 以 使 三 股 流体 的 总 质量 
流量 最 小 或 两 个 换 热 器 的 总 换 热 面积 最 小 . 
利用 通过 能 量 流 模 型 建立 起 的 换 热 系统 中 各 参 
数 间 的 约束 关系 ， 结 合 实际 设计 时 所 需 满 足 的 其 他 
约束 方程 ， 可 以 将 上 述 工程 优化 问题 转化 为 条 件 极 
值 问题 . 例如 , 以 三 股 流体 的 总 质量 流量 最 小 为 优化 
目标 , 结合 约束 方程 , 可 以 构建 以 下 拉 格 朗 日 函数 . 


F=mti+mati+mt+A(Aii+ A — 


Atotal) 十 


N N 
和 2 > 已 这 一 ou 十 入 3 = 2,i 一 Qun 十 


ji et 
Th 1—Tii— Qu (Rinit+ Rie1t+ Ri,al)— 

4 
Q2,1 (R2,a1t Ro,g,1+ R21)+e1,ate2,a 


N_1 
R 0 Rh 2 ai 
064 (Rihit+ Ri,giit Ri1,ad)— )+ 


这 了 Qui (Rupait Riginit Ri,aitl) 


bt Q2iR2,ait Rog,it R21 — (31) 
Q2an (Rezaitit Ro,ginlt R21in) 


其 中 , 入, wi 和 wp; 是 拉 格 朗 日 乘 子 . 
通过 拉 格 朗 日 函数 对 各 设计 变量 及 拉 格 朗 日 乘 

子 求 偏 导 ， 并 令 其 为 零 ， 可 以 获得 流体 的 总 质量 流 

量 最 小 时 ， 每 个 换 热 絮 的 最 佳 换 热 面积 和 各 文 路 流 

体 的 最 优 流量 所 需 满足 的 优化 方程 组 
OF/O0Q1i; =0,i=1,2,.…,N 
OF/0Q2i; =0,i=1,2,.… ,NN 
OF/Om:; = 0,i=h,a,l 
OF/0A; =0,i= 1,2 
OF/ON; = 0,i = 1,2,3,4 
OF/Ow; =0,i=1,2,... ,Nl1 
OF/0p; =0,i=1,2,... ,Nl1 


(32) 


上 述 优 化 方程 组 共 包 含 (4N +7) 个 方程 , (4N 十 
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7) 个 未 知 量 ,求解 该 方程 组 ， 可 以 直接 获得 最 优 的 
结构 和 运行 参数 , 从 而 实现 换 热 器 网 络 的 性 能 优化 。 
3 能 量 流 模型 的 分 析 及 优化 应 用 

3.1 相 变 换 热 器 的 性 能 分 析 

对 于 冷 流体 发 生 相 变 的 逆流 换 热 器 , 若 1) 高 
低温 工 质 的 进口 温度 、 质量 流量 和 比 定 压 热 容 , 2) 工 
质 相 变温 度 和 相 变 潜 热 ,3) 换 热 需 的 换 热 面积 和 对 
应 工 质 不 同 状态 时 的 换 热 系数 已 知 ， 通 过 联 立 求解 
方程 组 (11) 中 的 NN 个 方程 即 可 求 得 每 一 段 换 热 单 
元 的 换 热量 , 进而 实现 换 热 器 的 性 能 分 析 计 算 。 

在 求解 过 程 中 , 由 于 每 个 换 热 单元 的 换 热 系数 、 
工 质 比 定 压 热 容 以 及 换 热 融 的 热流 和 温度 分 布 存在 
耦合 关系 , 因此 需要 进行 迭代 计算 . 首先 , 设 定 每 个 
换 热 单 元 中 工 质 比 定 压 热 容 和 换 热 系数 的 初 值 ， 计 
算 每 个 换 热 单元 的 换 热 量 ， 根 据 求 得 的 热流 和 温度 
分 布 对 每 个 换 热 单元 中 工 质 的 比 定 压 热 容 和 换 热 系 
数 进行 修正 后 再 次 迭代 计算 . 重复 上 述 过 程 直到 每 
个 换 热 单元 中 工 质 的 比 定 压 热 容 和 换 热 系数 收敛 . 

对 换 热 单元 的 换 热 系数 和 工 质 比 定 压 热 容 进行 
修正 的 具体 方法 为 : 若 从 第 ;个 到 第 NN 个 换 热 单元 
的 总 换 热量 小 于 预 热 低温 流体 所 需要 的 热量 ， 则 更 
新 第 i 段 换 热 单元 内 低温 流体 的 比 定 压 热 容 为 工 质 
的 液态 比 定 压 热 容 ， 同时 该 换 热 单元 的 换 热 系数 更 
新 为 换 热 器 预 热 段 对 应 的 换 热 系数 数值 ; 若 该 换 热 
其 大 于 等 于 预 热 低温 流体 所 需要 的 热量 , 但 小 于 预 
热 低温 流体 需要 的 换 热量 和 低温 流体 发 生 相 变 需要 
的 潜 热 之 和 , 则 更 新 流体 的 比 定 压 热 容 为 无 穷 大 , 同 
时 换 热 系数 更 新 为 换 热 器 相 变 段 对 应 的 换 热 系数 数 
值 ; 若 该 换 热量 大 于 等 于 预 热 低温 流体 需要 的 热量 
和 低温 流体 发 生 相 变 需要 的 潜 热 之 和 ， 则 更 新 工 质 
的 比 定 压 热 容 为 气态 比 定 压 热 容 ， 同 时 换 热 系数 更 
新 为 换 热 器 过 热 段 对 应 的 换 热 系数 . 

对 于 换 热 面积 为 15 m? 的 道 流 相 变换 热 器 ,， 当 
高 温 工 质 的 质量 流量 为 50 kg/s， 入 口 温度 为 35 °C， 
比 定 压 热 容 为 1500 J/(kg:“C); 低温 工 质 的 质量 流 
若 为 10 kg/s， 入 口 温 度 为 15 "C， 相 变温 度 为 17 
"C, 相 变 潜 热 为 5000 J/kg, 液态 比 定 压 热 容 为 4200 
J/(kg. “0), 气态 比 定 压 热 容 为 2000 J/(kg:“C), 预 热 
段 、 两 相 段 、 过 热 段 的 换 热 系数 分 别 为 600、800、400 
W/n2…"C) 时 , 利用 本 文 构建 的 热 阻 网 络 模型 ， 取 
分 段 数 NN 为 100, 表 1 给 出 了 换 热 性 能 分 析 结 果 . 
将 计算 结果 分 别 代 入 换 热 器 预 热 段 、 相 变 段 和 过 热 
段 对 应 的 传 热 方程 和 能 量 守 恒 方 程 ， 各 个 方程 均 能 
得 到 较 精确 的 满足 ， 从 而 验证 了 本 文 提出 的 相 变 换 
热 器 热 阻 网 络 模型 的 准确 性 . 
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表 1 逆流 相 变 换 热 器 性 能 分 析 计 算 结 果 


Table 1 The performance solutions for the phase-change counter-fow heat exchanger 


枝 热 段 。 相 变 段 。 过 热 段 。 冷 流体 热流 体 

物理 量 “ 换 热 量 / 。 换 热 量 / 。 换 热 量 / 。 出口 温 ”出 口 温 
kJ kJ kJ 度 /*C 度 /*C 

84.41 49.94 22:77 18.15 32.90 


3.2 相 变 换 热 系统 的 性 能 优化 

对 于 图 4 所 示 的 换 热 系统 , 总 换 热 量 为 250 kW， 
总 换 热 面积 为 50 m?*; 高 温 工 质 的 入 口 温度 为 40 "C， 
比 定 压 热 容 为 5000 J/(kg*0), 低温 工 质 入 口 温 度 为 
4°C, 比 定 压 热 容 为 3000 J/(kg:“C); 中 间 回 路 工 质 的 
液态 比 定 压 热 容 为 3500 J/ (kg “C), 液态 密度 为 1000 
kg/m2, 气态 比 定 压 热 容 为 2000 J/ (kg“C); 储 液 镀 高 
度 为 2 m, 重力 加 速度 为 9.8 m/s?, 忽略 中 间 工 质 的 
状态 变化 对 换 热 系数 的 影响 , 为 1000 W/(m2:"C). 通 
过 合理 分 配 两 个 换 热 器 的 换 热 面积 和 各 回路 中 工 质 
的 质量 流量 可 以 使 三 条 回路 中 工 质 的 总 质量 流量 最 
小 . 求解 优化 方程 组 (32) 可 获得 该 换 热 器 网 络 在 上 
述 设计 条 件 下 的 最 优 结构 和 运行 参数 . 在 求解 过 程 
中 , 需要 解决 1) 中 间 回 路 工 质 的 相 变 温度 和 相 变 潜 
热 的 确定 问题 , 以 及 2) 中 间 回 路 工 质 比 定 压 热 容 和 
换 热 器 热流 分 布 的 耦合 问题 . 

中 间 回 路 工 质 的 相 变 温度 和 相 变 潜 热 取决 于 工 
质 在 两 换 热 器 中 的 压力 大 小 . 若 忽 略 中 间 回 路 工 质 
的 流动 阻力 , 则 泵 所 需要 提供 的 压 差 Ap 表示 为 

Ap = pgH — pgh (33) 

则 中 间 回 路 工 质 在 两 个 换 热 器 中 的 压力 p 可 以 由 下 
式 确定 : 


p=palt+pgh+i+ Ap= pait+pgH (34) 
其 中 ,pai 为 储 液 负 中 液体 温度 所 对 应 的 饱和 蒸汽 
压 ,， 互 为 储 液 缸 高 度 , /为 储 液 负 中 液 面 高 度 . 由 
以 上 分 析 可 见 ， 中 间 回 路 工 质 的 压力 受 其 在 换 热 器 
2 中 的 出 口 温度 影响 。 同 时， 压力 的 大 小 决定 了 工 
质 在 两 个 换 热 器 中 的 相 变 温度 和 相 变 潜 热 ， 进 一 步 
会 影响 到 各 换 热 器 内 工 质 的 比 定 压 热 容 和 热流 分 布 ， 
继而 又 影响 到 换 热 器 2 的 出 口 温 度 . 因此 , 在 该 闭 
式 换 热 系统 中 ， 工 质 的 工作 压力 和 换 热 存在 着 复杂 
的 看 合影 响 关 系 , 需要 进行 近代 求解 . 

表 2 给 出 了 相 变 换 热 系统 的 优化 结果 ， 包 
括 每 条 支 路 流体 的 质量 流量 和 每 个 换 热 器 的 面积 
大 小 . 

图 6 给 出 了 优化 后 各 支 路 工 质 流量 和 换 热 句 面 
积分 配 随 换 热 量 的 变化 关系 . 当 总 换 热 面积 一 定时 ， 


次 流体 。 准 流 体 
预 热 段 。 相 变 段 过 热 段 。 相 变 起 。 相 变 终 
物理 量 ” 换 热 面 换 热 面 换 热 面 始点 热 。 止 点 热 
积 /m2 积 /m? 。 积 /m2 。 流体 温 流体 浊 
度 /C_ 。 度 /c 

计算 
结果 8.10 3.60 3.30 34.03 34.69 


各 文 路 工 质 流 其 随 换 热量 的 增 大 而 提高 ， 且 提高 速 
度 变 快 ; 换 热 量 越 大 , 换 热 面 积分 配 越 均 匀 . 


表 2 相 变 换 热 系统 的 性 能 优化 结果 
Table 2 The optimal results for the heat 


exchange system involving fluid phase changes 


mn/kg:s™!l ma/kgs-! m/kgs™! Ai/m? A2/m? 
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图 6 换 热 系 统 优化 设计 参数 随 换 热 负荷 的 变化 图 


Fig. 6 The optimal mass fow rates of working fluids and the 
allocations of thermal conductance versus the heat 


transfer rate 


8 
5 ea a 
et 
4 ~ oo J0.96 
oo 
es ae。 | 
3 口 
有 0.95 
40 42 44 46 48 30 52 54 56 58 60 
Aiva/m? 


图 7 换 热 系统 优化 设计 参数 随 总 换 热 面积 的 变化 图 
Fig. 7 The optimal mass flow rates of working fluids and the 
allocations of thermal conductance versus the total thermal 
conductance 
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图 7 给 出 各 优化 参数 随 换 热 器 总 换 热 面积 的 变 
化 关系 . 当 换 热 负 荷 一 定时 ， 各 文 路 工 质 流量 随 总 
换 热 面积 的 增 大 而 减 小 ， 而 总 换 热 面积 越 大 ， 换 热 
器 的 换 热 面积 优化 分 配 越 不 均匀 


4 


工 质 在 换 热 过 程 中 发 生 相 变 时 , 对 工 质 物 性 、 换 
热带 换 热 系数 和 闭 式 换 热 系 统 的 工作 压力 均 产生 影 
响 ， 加 剧 了 系统 性 能 分 析 及 优化 的 复杂 度 . 本 文 在 
考虑 工 质 发 生 相 变 的 前 提 下 ， 建 立 了 顺 流 、 逆 流 和 
又 党 三 种 典型 换 热 器 的 等 效 热 阻 网 络 图 (能 基 流 模 
型 )， 获 得 了 换 热 器 性 能 与 其 结构 /运行 参数 和 工 质 
物性 间 的 直接 物理 联系 ， 并 通过 对 逆流 相 变换 热 需 
的 性 能 计算 分 析 ， 验 证 了 本 文 提出 的 热 阻 网 络 模型 
的 准确 性 . 

以 此 为 基础 , 针对 含 工 质 相 变 过 程 的 换 热 系统 ， 
构建 了 反映 系统 中 热量 整体 输 运 规律 的 能 量 流 模型 ， 
结合 基 尔 霍 夫 电压 和 电流 定律 建立 了 各 独立 设计 变 
量 间 的 整体 约束 关系 ， 并 以 多 回路 相 变换 热 系 统 为 
例 ， 应 用 拉 格 明日 乘 子 法 获得 了 工 质 总 质量 流量 最 
小 时 的 系统 最 优 结 构 和 运行 参数 所 需 满足 的 控制 方 
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程 组 . 求解 该 控制 方程 组 ， 实 现 了 相 变 换 热 系统 换 
热 性 能 的 整体 优化 . 
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